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O . Abstract 

.^^ A new idea for turbulent skin-friction reduction is proposed, wherein 

the shape of the sohd wall is designed to create the spanwise pressure 
>D gradient acting similarly to the well-known method of drag reduction by 

' in-plane spanwise wall motion. Estimates based on the assumption of 

^.^ similarity with drag reduction effect of in-plane wall motion suggest that 

r"| drag reduction of about 2.4% can be achieved in the flow past a wavy wall, 

^-^ with the crests forming an angle of about 38° with the main flow direction, 

■ O and the wave period in the main flow direction equal to about 1500 wall 

units. The required height of the wall waves have to be adjusted to achieve 

the same intensity of the spanwise motion as that created by an in-plane 

moving wall of the same wavelength and with peak wall velocity equal to 

^ 2 wall units. Further research is being conducted in order to determine 

. ^ this height. Suggestions are made for further research on confirming the 

C^ feasibility of the proposed method and on optimising the wall shape. 

, ^ 1 Introduction 



I 



Turbulent skin-friction control techniques are classified into passive and active. 

^ Passive techniques do not require energy input. For example, the shape of the 

soUd surface can be made such that the skin friction will be less in the flow past 

this surface than the skin friction in the flow past a flat wall at the same condi- 
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^^ tions. The only surface shape known, reliably, to decrease drag is one covered 



with riblets '11 . Riblets are longitudinal grooves in the surface exemplified in 
/--s Figure [T^. Riblets inhibit spanwise velocity fluctuations, that are fluctuations 

(T^ directed across the grooves, and this modification of the structure of turbulence 

reduces the drag. The cross-sectional shape of the riblets can be different from 
that shown. Importantly, to reduce the drag, both the distance between the 
neighbouring groves and their height should be about 15 wall units. A wall unit 
is defined as the distance normalised by the wall skin friction and fluid viscosity. 
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S For aircraft applications a wall unit might be of order of one micron, and ri- 

blets have, therefore, to be very small. The drag-reduction level achievable with 
riblets is less than 12%, and the extra manufacturing and maintenance costs 
involved make riblets only marginally effective in a practical environment. 

Active control techniques require energy input. For example, the surface of 
the wall can perform oscillatory in-plane motion in such a way that the turbu- 
lent skin-friction drag is reduced. Spanwise oscillations of the wall can reduce 
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Figure 1: a) Riblets; b) Spanvise wall motion. 



the skin- friction drag by up to 40% [2j . This effect has been predicted computa- 
tionally and confirmed experimentally [3]. More complicated in-plane motions, 
in which different parts of the wall surface move with different time-dependent 
velocities, can produce even higher effects as was shown in 0]. Importantly, 
the energy required for such an actuation need only be about half of the gain 
obtained due to the drag reduction, so that the energy budget is favourable. Of 
particular interest for the present work are the findings of ^5,, where a steady 
wall motion, illustrated in Figure[T]D, was considered. The authors showed that 
net energy savings of 23% were possible. They also found that the optimal 
longitudinal wavelength of the forcing is somewhat larger than 1000 wall units. 
Other approaches to turbulent-flow control, such as blowing and suction, us- 
ing micro-electro-mechanical actuators with feedback control, or using plasma 
actuators for creating the cross-flow motion, have been proposed, but so far 
none of these proposals has resulted in practical applications. With the above 
described approaches not being widely used in practice, there remains a need 
for a practical and cost-effective control system allowing large skin-friction re- 
ductions to be achieved. The present work proposes a simple and practical 
method of passive control intended to achieve the same effect as the spanwise 
wall motion shown in Fig. [lb. 

2 The idea 

The idea consists in generating, by selecting the appropriate shape of the surface, 
the cross- flow motion producing the drag- reduction effect. An example of such 
a wall shape is shown in Figure[2J The deflection of the main flow by the wall 
creates different spanwise pressure gradients on the upwind and downwind sides 
of the smooth waves. This pressure gradient pushes the fluid sideways, similar 
to a spanwise wall velocity shown in Figurellb. 

Unlike the case of active control, the rigid wall will not require energy input. 
However, the additional motion generated by the non-fiat wall will result in an 
increase in energy dissipation similar to that for a moving wall. This additional 
energy dissipation will manifest itself as a pressure drag on the non-flat surface. 
Therefore, the proposed method can be considered as a simple and feasible ap- 
proach to providing energy for generating the drag-reducing spanwise motions. 

Crucially, even a relatively small variation of the wall shape will produce 




Figure 2: Drag-reducing wavy wall 



a significant variation in the velocity near the wall, where the drag-reduction 
mechanism is concentrated, because of the mechanism involved in the well- 
known phenomenon of viscous- inviscid interaction (see, for example, [6]). The 
properties of boundary layer flow are such that the displacement of streamlines 
caused by variation in wall surface shape is passed with little change to the 
streamlines at the outer edge of the boundary layer, where the flow velocity is 
large. This displacement then results in the pressure variation proportional to 
the velocity at the outer edge of the boundary layer and to the displacement 
magnitude. The pressure variation is then passed back to the wall, since the 
pressure does not vary significantly across the boundary layer. Near the wall 
this pressure variation results in significantly larger variation of the velocity, 
because the velocity itself is small near the wall. At the intuitive level this can 
be understood by recollecting the well-known Bernoulli equation, from which 
it follows that the pressure variation and the variation of velocity squared are 
of the same order of magnitude. At the outer edge of the boundary layer the 
velocity variation At/, created by a small displacement of the streamlines, is 
small as compared to the velocity U itself, so that the pressure variation Ap is 
of order UAU, but near the wall, where the velocity approaches zero, the velocity 
variation Au is of the same order as the velocity itself, so that the same Ap 
corresponds to Au^. Hence, Aw is of the same order as {U/Au)AU, that is much 
greater than AU. For this effect to take place, the characteristic thickness of the 
near-wall layer should be small as compared to the longitudinal and spanwise 
length scale of the wall-shape variation. Fortunately, this is possible since the 
characteristic scale of the near-wall layer is about 100 wall units, as evidenced 
by the typical streak spacing, while the optimal longitudinal wavelength of the 
steady spanwise wall motion (FigurefTl) is somewhat greater than 1000 wall units 
— that is, one order of magnitude greater. This wavelength is almost two orders 
of magnitude greater than the wavelength of riblets. Hence, it is much easier to 
manufacture and maintain in practice. 

Direct numerical simulation of a flow past such a wall is more computation- 
ally expensive than a simulation of a flow past a flat wall, not only because of a 
more complicated geometry, but also because the waviness of the wall adds an 
additional spanwise length scale, which is about an order of magnitude greater 
than the streak spacing. Before performing expensive direct numerical simu- 
lations one can attempt to use the similarity between the action of spanwise 
pressure gradient and the action of spanwise wall motion in order to determine 
the range of the parameters where the drag reduction is most likely to be ob- 
served. This similarity was already reported in [5] where drag reduction by 
wall oscillations was first demonstrated. In the case considered in that paper 
the similarity was in fact exact, since the two situations could be made iden- 
tical by a simple change of the frame of reference. In our case this can only 



be approximate. We will assume that if the spanwise motion generated by a 
particular wavy wall is close to the spanwise motion generated by a particular 
in-plane spanwise wall motion then the effect of this two methods of control 
will be similar, resulting in the same reduction in the skin friction at the wall. 
The difference between the two methods of control is that the spanwise motion 
of the wall consumes power, while a rigid wavy wall does not consume power. 
However, the wavy wall will experience pressure drag, not present in the flow 
past a spanwise moving wall, and overtaking this additional drag will require 
additional power. Thus, in both cases there is a price to be paid for control, 
which can be expressed as the power required for control. Viotti, Quadrio, and 
Luchini [S] calculated both the power saved, Pgav and the (negative) power re- 
quired, Preq, for the spanwise-moving wall, and considered their difference, Pnet, 
which characterised the net gain achievable by the method in principle. We will 
find the shapes of the wavy walls generating a spanwise shear similar to the 
spanwise shear obtained in several cases by |5], and assume that Psav is the 
same in both cases. We will then calculate Preq for a flow past a wavy wall, us- 
ing an approach demonstrated in |5] to work well for the spanwise-moving wall, 
thus obtaining Pnet for the wavy wall. We then optimise it to get an estimate 
of the drag reduction achievable in the flow past a wavy wall and the optimal 
parameters of the wavy wall. 

3 Boundary layer on a wavy wall and the power 
required 

3.1 Boundary layer equations for the wavy- wall case and 
the spatial Stokes layer case 

Similar to |5], the boundary layer equations, linearized around a linear proflle 
will be used. Written in wall units, they have the form 
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3.1.1 The spatial Stokes layer case 

In the spatial-Stokes-layer (SSL) case the flow is driven by the movement of the 
boundary, so that u — v — Q, w — We^^^^ at y = 0. The pressure gradient is 
zero, and Uss\ = Vssi = 0. Taking w — Wwss\{y)&^'^''^ ^ where y — kx y gives 

(Pwss\ ,,. 

tywssi = ^-2 (1) 



with boundary conditions lissi = 1 at y = 0, and Wssi — > as y — > oo, the same 
as in fS! . We solved (IT]) numerically and obtained a perfect agreement with the 
formula Wssiiy) = Ai{—iye^'^^^^^)/Ai{0), equivalent to (7) in [5]- 

3.1.2 Wavy wall case 

In the wavy wall case the boundary condition isu = v = w = Oaty = 0. 
The pressure distribution is generated by the inviscid flow above the boundary 
layer. We assume it to have the form p — pe'*^''="^+'^^^'. Accordingly, {u,v,w) — 
(M(y),'0(y), w(y))e'(''"^+''"^). This leads to 

ikxyii + V — —ikxP + u" (2) 

ikxyw = —ikzP + w" 
ikxii + v' + ikzW = 
Eliminating v and then rescaling as 

w{y) = ■ikzk~'^^^pw{y), u{y) = iklk^^/'^pu{y) 

gives, after simple transformations, the following system 
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ly— -iw= —^ 
dy dyi 

d'^w 

with boundary conditions ^ — > and w — )• as y — > oo, and u = w = 0, 
it" = k'^/k'l at y = 0. The last condition is in fact ([2| taken at the wall. The 
equation for w separates from the equation for u and can be solved first. It is 
then natural to take 
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with u^ and itp satisfying 



.duyj . _ d^u„ 
dy dy^ 



ly— — iw — ,,„ , Uu, — >■ as y — >■ oo, and Uw = U, u.,, = 1 at y = 0, 



iy—r3- — ,,i' , Up -^ as y -> oo, and u„ = 0, u" = at y = 0. 
dy dy-^ ^ f 

From the physical viewpoint ityj corresponds to the perturbation of u due to 
wall- normal velocity induced by spanwise velocity dependence on z, while Up 
is related to the perturbation of u due to the longitudinal pressure gradient 
induced by the wall. 

These ordinary differential equations were solved numerically using Mathe- 
matica. 







Figure 3: Comparison of spanwise shear RealWgj,[(y)e*'^^'^ (solid) with 
RealC™u;'(2;)e*('=-+'^") (dashed), for Kx/{2tt) = 0,1/6,2/6,3/6,4/6, and 5/6. 



3.1.3 Matching the spanwise shear profiles in the SSL and wavy- wall 
case 

At the wall the spanwise velocity is zero in the wavy-wall case and nonzero (and 
has a maximum amplitude) in the SSL case. Therefore, they cannot be directly 
matched. However, a simple translational motion might be not so important 
because of Galilean invariance. We presume that it is the spanwise shear that 
favourably affects turbulence leading to drag reduction. Therefore, we seek 
to match the SSL-case spanwise shear Ww'^^^{y)e''^'''^ with the wavy-wall-case 
spanwise shear ik^kx pw'(y)e*'^'^="^"'"*'^^^. Note that the wavy-wall period in 
the spanwise direction is expected to be much larger than the characteristic 
scale of near-wall turbulent structure and, hence, dependence on z can simply 
be neglected. On the other hand, a phase shift between the SSL and wavy- 
wall cases is acceptable. Accordingly, the matching was done by minimising 
numerically 
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over C = ikzkx p/W and (p. This gave C = C„ 



0.8980 and 



1.5708. Hence, to achieve matching, the height of the wall waves should be such 

that 

p = C„,iy/(zfc,fc-2/3) (3) 

Note that 0„i w tt/2. Figure Is] shows the quality of the matching. 

We believe that the agreement is sufficiently close to expect that these two 
spanwise profiles will lead to similar magnitudes of skin-friction reduction. For 
completeness, the comparison of the spanwise velocity profiles is shown in Fig- 
ure ll 



3.1.4 Comparison of the power required for control in the SSL and 
wavy-wall cases 

The power per unit wall area (on one wall) required to drive the SSL flow can 
be expressed in wall units as 
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Figure 4: Comparison of spanwise velocity Realwssi(i')e''''^'^ (solid) with 
RealC™u;(y)e^('=-+'^'") (dashed), for Kx/{2tt) = 0, 1/6,2/6,3/6,4/6, and 5/6. 



while the corresponding expression in the wavy-wall case is $+ — $^J; -|- <f>^ 
where 

/"OO 

^t,=k,k-^'^p |zD'(y)|V2dy, 



Kk-^'^ / \u'{v)\^/2dv. 
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lip is selected by the matching rule (|3| then <&+ < $j^[. This is an artefact 
of the specific way matching was done: optimising over C is equivalent to pro- 
jecting the SSL solution onto the direction of the wavy-wall solution in the L2 
functional space: the length of the projection of a vector is always less than or 
equal to the length of the vector itself. To be on the safe side, therefore, it is 
better to assume that the energy needed to drive the spanwise component of the 
wavy-wall flow to achieve the same skin friction reduction as in the SSL flow is 
the same as the energy needed to drive the SSL flow. With this assumption, if 
the equivalent wavy-wall flow had no u component it would generate the same 
overall net drag reduction as the SSL flow. In reality, the u component, however, 
results in an additional energy dissipation, $+. Hence, the power required to 
generate the wavy-wall flow equivalent to the SSL flow will always be greater, 
with the ratio equal to 



5i 



\w 



1^2 Up 



w^W 



1,2 



7,1 12 



11*11 



Here, the squared norm ||.|p = J^ |.p dy. The ratio r turns out to be dependent 
on k,j./kz only, and it is also obvious that there exist a value of k^/k^ for which 
the power required will have a mininuim. Calculating the norms numerically 
gives 

t2 l2 

r = 3.122 + 2.323-f + 0.7986-^. 
Minimizing r gives r^in = 5.846, which is attained at 



= 0.7657. 



opt 



This corresponds to the angle between the mean flow direction and the 
direction perpendicular to the wall crests and troughs 0opt = 52.56°. 
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Figure 5: Power required to move the wall in SSL flow as a function of the 
wavelength A+; W+ = 6 (soHd), W+ = 12 (dashed), DNS [5, (dots). 



3.2 Drag reduction estimate 

In [5] the drag reduction was characterised by the net gain in the power needed 
to drive the flow, expressed as a percentage share of the power needed to drive 
the uncontrolled flow. The power needed to drive the uncontrolled flow in a 
plane channel, if expressed in wall units, is 
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where Uj^ is the bulk velocity. We will use the same empirical formula for the 
skin friction coefficient Cj = 0.0336i?e^"'^ ^ as in [5 . The net power gain P^ct 
was obtained as a sum of the power Pgav saved due to the reduction of the 
skin friction and the (negative) power Pjcq required to drive the in-plane wall 
motion. In [5 Pgav and P^cq were obtained from direct numerical simulations, 
and expressed in % of the power needed to drive the uncontrolled flow, and the 
same convention is used here. It was also shown that with good accuracy Prcq 
can be also obtained from the solution of the linearized boundary layer equation 
(fTl). It should be understood, of course, that (IT]) is formulated in wall units of 
the controlled flow while Pj-cq is presented in 5J as a function of the wavelength 
Aj in the units normalised with the skin friction of the uncontrolled flow. The 
reduction of the skin friction is proportional to 100% — Psav, which, after some 
manipulation, gives the formula 
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This formula, although not given explicitly in [5], was obviously used there. 
According to our calculations and the solution given in [S], J^ |Wgj,[(y)p/2dy = 
0.3157. Comparison of this formula with the DNS results of [S], measured from 
a digitised plot in that paper, is shown in Figure [5] and is in fact equivalent to a 
part of their Figure 6. The agreement confirms that we use the same approach 
as in [S]. 

Using the data provided by the authors of [S] , the following fits were obtained 
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Figure 6: Power saved fits (ill? I, bottom to top, for the DNS data of [5]. 



for Psav for W+ = 1,2,6 and 12: 

P,^^ 1 = 1.135 + 0.002929A+ - 1.205 • 10"^A+^ + 1.447 • 10"^°A+^ 

- 1.047 • 10~"A+^ + 2.609 • 10"l^A+^ (4) 



Psav,2 = -1.856 + 0.03954A+ - 5.28537 • 10"^A+^ + 3.498 • lO^^A 

- 1.127 • 10""A+^ + 1.328 • lO^^^A 



5 

X 7 



(5) 



P^^^Q = 15.25 + 0.04888A+ - 4.441 • 10"^A+^ + 1.628 • 10~^A+^ 

- 2.845 • 10^12 A+^ + 1.938 • lO^^^A 
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X 7 



(6) 



P 



sav,12 



27.90 + 0.03824A+ - 2.810 • 10"^A^ 



1.082 •lO"^^ A: 



;.015-10"9A+ 
:V 5.535 -lO^^^AJ 



(7) 



The quahty of the fits is illustrated by Figure [6j 

When the wavy wall is used to create the spanwise shear instead of the 
in-plane wall motion, the power required should be multiplied by the ratio r 
depending on the angle of the wall wave, with the minimal value r^- The net 
power saving is then Pgav + ?'Prcq- Figure [T] shows the result for r — r,„ for two 
wall heights equivalent to W^ = 2 and 6, together with the W^ — 6 SSL case. 
As one can see, in the case corresponding to W^ = 2 the wavy wall can be 
expected to give a positive net power saving. The maximum of about 2.4% drag 
reduction is attained at A+ « 1520. In the case corresponding to a wavy wall 
equivalent to SSL with W^ — 6 only drag increase is predicted. We interpolated 
between W^ — 1, W'^ = 2, and W'^ = 6 cases, and it turned out that W^ = 2 
gives very nearly the best result. 

Note that the calculations were done using the Mathematica package, and 
while the parameters affecting the accuracy were varied to verify it, the actual 
accuracy might be somewhat less than four digits given here. In any case, the 
nature of the present study is such that the values obtained are only indicative. 
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Figure 7: Net power savings for wavy wall case: Psav,2 + r„i PYcq\]y+^2 (solid 
line) and Psav,6+^m -Proqlvy+^g (thick dashed line) and for SSL: DNS for W^ = 6 
(dots) and Psav,6 + -Prcql^y+^g (thin dashed line). 

4 Discussion and conclusions 

The present analysis is based on assumptions which, although being reasonable, 
have not yet been confirmed. The most natural next step will consist in per- 
forming direct numerical simulation of a flow past a wavy wall and checking that 
the drag is indeed reduced. This task is more complicated than a now-standard 
direct numerical simulation of a flow in a flat channel, not only because of the 
more complicated geometry, but also because the wavy wall introduces an ad- 
ditional length scale in the spanwise direction, which is noticeably greater than 
the typical spanwise scale of near-wall structures. For this reason the spanwise 
dimension of the computational domain might need to be increased. The effect 
itself is not particularly large, and it is expected to be observed in a relatively 
narrow range of wall shape parameters. Therefore, a trial-and-error approach 
to selecting these parameters in direct numerical simulation could be inefficient. 
Providing a reasonable estimate for these parameters was the main motivation 
of the present study. 

Note that this study is based on the direct numerical simulation results [S] 
for Rct- = 200, and the conclusions for other values of the Reynolds number 
might be somewhat different. 

Concerning the required wave height, our results only indicate the magnitude 
of the spanwise motion that the wall should generate. This might be enough 
for direct numerical simulations, since estimating approximately this magnitude 
for a particular wall height does not require averaging over a long time: here, 
trial-end-error is more likely to work. From the mathematical viewpoint, the 
required estimate could be obtained by finding the inviscid outer limit of the 
asymptotic expansion of the flow past a wavy wall as the Reynolds number tends 
to infinity. This was not undertaken because first direct numerical simulations 
are likely to be tried for relatively small Reynolds numbers. An alternative is 
to perform an analogue of the analysis presented in this work numerically using 
an eddy viscosity model. This line of research is being pursued currently by the 
colleagues of the author. 

Drag reduction by 2.4% is much smaller than what can be obtained by 
in-plane wall motion. However, the method proposed here is much easier to 
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implement in practice, since it is passive and does not require wall motion. It is 
also easier to implement than riblets, because the wavelength of the proposed 
wall is two orders of magnitude greater than the wavelength of riblets, and 
because the height-to- wavelength ratio for the proposed wall is also much smaller 
than that of riblets. 

The present analysis is based on the assumption that similar distributions 
of spanwise shear will result in similar decrease in turbulent friction, while the 
longitudinal perturbation of the mean velocity will not affect turbulent friction. 
This can only be verified in direct numerical simulations or experiment. Note 
that for the very first tests numerical simulations need not be done for a wavy 
wall; the same effect might be expected to be achieved if a steady body force is 
applied simulating the effect of the pressure distribution caused by a wavy wall. 

Another significant source of inaccuracy in the above analysis is the use of 
linearized boundary layer equations and the assumption that the mean profile 
in the uncontrolled flow is linear. Therefore, the particular value of 2.4% should 
be considered as indication of the accuracy with which the drag needs to be 
determined in direct numerical simulations rather than a definite prediction. 
Similarly, the optimal parameters of the wall obtained here are just an indication 
of the reasonable initial approximation for the search of the true optimal. 

It might be possible to increase the effect by optimising the wall shape, since 
it need not be of the form sin {k^x + k^z). More complicated shapes might be 
better. Optimisation over complex wall shapes would be impossible on the ba- 
sis of direct numerical simulations with computing power currently available. 
Fortunately, recent studies [71 H] suggest a much more efficient than direct nu- 
merical simulation, albeit approximate, approach to such optimisation. Once 
approximate optimisation using these approaches has been performed, further 
refinement might be done by direct numerical simulations and experiment. 

Since the riblet wavelength and the proposed wall wavelength are so different, 
they can be combined to work simultaneously. It has already been demonstrated 
that the effect of wall oscillations and effect of riblets are almost additive [9] . 

The main conclusion of the present work can be formulated in the following 
way. There are reasons to expect that a turbulent flow past a wall the height 
of which is of the form h — H sm{k^x^ + k^z), where x^ is the coordinate in 
the main flow direction, with fc+ = 27r/A+ « 27r/1520 and k^ /k^ « 0.7657 and 
with a suitable H can be expected to have the skin friction about 2.4% less than 
the flow past a flat wall, other things being the same. The value of H should 
be such as to generate spanwise velocity of the order of 2 wall units. The wall 
units here are based on the skin friction in the flow past a flat wall. The values 
given are only indicative. 

This work would be unlikely to appear without the research environment 
provided by the large team working on the project on turbulent drag reduction 
under the grant EP/G060215/1, funded by EPSRC together with Airbus Oper- 
ations Ltd and EADS UK Ltd. The author would like to use this opportunity 
to thank all his colleagues in this team. 
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